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This  report  is  an  attempt  to  analyze  an  engagement  be- 
tween an  aircraft  and  a ground  threat,  each  of  which  is 
using  optics  in  tracking  or  electronic  countermeasures 
roles.  Actually,  the  report  is  limited  to  staurting  the  an- 
alysis, which  will  probably  continue  far  into  the  future. 

In  limiting  the  scope  of  the  problem,  the  atmospheric  ef- 
fects have  been  ignored.  Also,  the  aircraft  has  not  been 
allowed  to  manuever  but  just  flies  in  a straight  path  over 
the  threat.  This  is  not  realistic,  but  necessary  in  order 
to  concentrate  on  the  real  problems  such  as  tracking  and 
jamming. 

Dusting  off  and  reviving  my  expertise  in  FORTRAN  pro- 
gramming was  a herculean  task  cheerfully  accepted  by  Mr. 
William  McQuay  of  the  Avionics  Laboratory.  I would  also 
like  to  thank  Capt.  Phil  Gordin  for  his  assistance  in  estab- 
lishing the  direction  and  scope  of  this  problem. 

Special  thanks  goes  to  Dr.  Shankland  for  his  guidance 
and  direction  in  pursuing  the  problem. 

I wish  to  acknowledge  my  indebtedness  to  my  wife  for 
keeping  things  under  control  at  home  and  helping  in  the  pre- 
paration of  this  report. 


James  R.  Penick 
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List  of  Definitions 


User  Variables 


ANGI£ 

ASPEDX.ASPEDY.ASPEDZ  . 

BARK 

BITE 



DELAY  

DOGl 


Used  in  FB4.  The  solution  to  the 
quadratic  formula  used  to  deter- 
mine the  point  of  intercept. 

The  speed  of  the  aircraft  in  the 
X,  y,  and  z directions.  Used  in 
all  Fly  Bys. 

Used  in  FB4.  Represents  a small 
part  of  the  a and  c constants  in 
the  quadratic  formula  used  to  find 
the  point  of  intercept. 

Used  in  FB4.  Another  step  up  from 
BARK  in  determining  the  a and  c 
constants. 

Used  in  FB4,  Represents  a change 
in  range.  Output  of  RNORM. 

Used  in  FB4.  Represents  the 
amount  of  time  needed  to  fire 
another  missile  after  one  has 
been  fired. 

Used  in  FB4.  Represents  the  a 
constant  in  the  quadratic  formula. 


D0G2 Used  in  FB4.  Represents  the  c 

constant  in  the  quadratic  formula. 

2 

DOG 3 Used  in  FB4.  Represents  b in  the 

quadratic  formula. 


Dog4 Used  in  FB4.  Represents  b in  the 

quadratic  formula. 

g Used  in  FB4.  Represents  a small 

change  in  elevation.  Output  of 
RNORM. 

EFFRM  .........  Used  in  FB4.  This  is  the  effective 

range  of  the  missile. 


EFPRT  This  is  the  effective  range  of  the 

tracking  radar.  Used  in  FB4. 

p Represents  a small  change  in  azimuth 

and  is  used  in  FB4. 


V 


HIALT 


Variable  used  in  all  Ply  Bys  to 
simulate  detection  probability. 
Divided  by  aircraft  range  to  ar- 
rive at  detection  probability 
number. 

JFLAG Used  in  FB4.  This  flag  is  set  in 

SCANT  and  used  to  determine  if  the 
tracker  is  detecting  the  aircraft 
for  the  first  time.  Statistics 
are  collected  only  when  the  tracker 
sees  the  aircraft  for  the  first 
time . 


KFLAG Used  in  FB4.  Used  in  DETECT  to 

determine  if  detection  by  the 
acquisition  radar  has  been  con- 
tinuous. 

MFLAG Used  in  FB4.  Used  in  DETECT  to 

determine  if  the  aircraft  is  being 
detected  by  the  acquisition  radar 
for  the  first  time. 

MINEPR This  is  the  minimum  effective  range 

of  the  missile.  Used  in  FB4. 

MFLTT Used  in  FB4.  This  is  the  missile 

flight  time  from  launch  to  inter- 
cept. 

NACFT Used  in  FB3  FB4.  This  is  the 

nvimber  of  aircr'ift  remaining  in 
the  simulation.  Used  as  a counter. 
When  0,  the  simulation  ends. 

NACFTI Used  to  record  the  number  of  air- 

craft at  the  beginning  of  the  sim- 
ulation. Used  in  FB3  and  FB4. 

NMISL Used  in  FB4,  This  is  the  number 

of  missiles  the  threat  can  fire 
at  a single  aircraft. 

NNSURR Used  in  FB4.  This  is  the  number  of 

aircraft  shot  dovm. 

NSURR  Represents  the  number  of  aircraft 

that  survive.  Used  in  FB4. 

PHE Used  in  FB4.  This  is  the  angle  of 

azimuth  used  by  the  tracker  to  com- 
pute the  aircraft  position  and  an 
intercept  point. 


Vi 


POSAX.POSAY.POSAZ 


Used  in  all  Fly  Bys,  Represents 
the  aircraft  position  at  TNOW, 


POSAXI , POSAYI , POSAZI  . 
POSMX.POSMY.POSMZ  . . . 

POSTX.POSTY.POSTZ  . . . 
PRBDET  

PRBTRK  

RVEL 


Used  in 
initial 


all  Fly  Bys.  This  is  the 
position  of  the  aircraft. 


Used  in  FB4.  This  is  the  position 
of  the  missile  at  the  time  of  mis- 
sile-aircraft intercept. 

Used  in  all  Fly  Bys.  This  is  the 
position  of  the  threat. 


Used  in  FB2,  FB3.  and  FB4.  Repre- 
sents the  probability  of  detecting 
the  aircraft  with  the  acquisition 
radar. 

Used  in  FB4,  This  is  the  proba- 
bility of  detecting  the  aircraft 
with  the  tracking  radar. 

Used  in  FB4.  This  is  the  speed 
of  the  aircraft. 


SCANR This  is  the  scan  rate  of  the  acqui- 

sition radar.  Used  in  FB2,  FB3t 
and  F34. 

SCANRT Used  in  FB4.  This  is  the  scan 

rate  of  the  tracking  radar. 

SS(1) Used  in  all  Fly  Bys.  This  is  the 

only  state  variable,  the  aircraft 
range. 

SST  The  aircraft  range  as  determined 

by  the  tracker.  Used  in  FB4. 

TANGLE This  is  the  masking  angle  of  the 


horizon.  The  acquisition  radar 
cannot  "see"  below  this  a.igle  due 
to  terrain  masking.  It  is  measured 
with  respect  to  the  threat  position. 
Used  in  FB3  and  FB4. 


TDETI Used  in  FB4.  This  is  the  time  of 

initial  detection  by  the  acquisi- 
tion radar.  It  is  reset  if  the 
detection  is  not  continuous. 

TEFFR Used  in  all  Fly  Eys.  This  is  the 

effective  range  of  the  acquisition 
radar. 
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THETA This  is  the  angle  of  elevation 

computed  by  the  tracker.  Used  in 
PB4. 

TTIME This  represents  the  total  time  of 

continuous  detection  by  the  acqui- 
sition radar.  Used  in  FB4. 

X A random  number  generated  by  BRAND. 

It  is  compared  to  PRBDET  to  deter- 
mine if  detection  has  occurred. 

Used  in  FB2,  FB3,  and  FB4. 

XMISS Used  in  FB4.  This  is  the  miss  dis- 

t£mce  of  the  missile. 

XMSPED This  is  the  missile  speed.  Used 

in  FB4. 

Y A random  member  generated  by  RNORM. 

Used  in  FB4  to  determine  if  de- 
tection has  been  made  by  the  track- 
ing radar. 


GASP  Variables 


ATRIB(X) This  is  a buffer  for  attribute 

values  being  stored  in  or  re- 
trieved from  QSET,  the  file  stor- 
age area.  Used  in  all  Fly  Bys. 

LFLAG(X) Used  in  all  Fly  Bys.  These  are  the 

state  condition  flags. 

MSTOP Used  in  all  Fly  Bys.  This  is  GASP 

indicator  which  specifies  how  the 
simulation  will  be  ended.  If 


MSTOP  is  read  in  as  a 0,  then  it 
must  be  reset  to  a negative  value  to 
stop  the  simulation.  If  read  in 
as  a positive  value,  then  the  sim- 
ulation will  continue  until  TNOW 
reaches  the  user  specified  end 
time , TTFIN . 

TNOW  .........  Used  in  all  Fly  Bys.  This  is  the 

current  time  of  the  simulation. 
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Abstract 

This  report  contains  data  on  the  construction  of  a com- 
puter model  to  study  a one-on-one  engagement.  It  is  built 
so  that  parameters  in  the  RF  or  optical  regime  can  be  used. 

The  aircraft  flies  over  the  threat,  is  not  allowed  to  maneuver, 
and  the  atmosphere  has  been  ignored. 

Fly  By  1 introduces  xhe  techniques  employed  by  the 
GASP  IV  simulation  language.  The  aircraft  is  detected  when 
it  comes  within  range.  In  Fly  By  2,  a more  probabilistic 
determination  of  detection  is  used,  and  the  radar  scans  for 
the  aircraft.  Fly  By  3 makes  20  runs  of  the  Fly  By  2 program 
and  collects  statistics  on  the  rauTge  of  detection.  Ply  By 
4 incorporates  a track  and  fire  role  into  the  threat.  It 
shoots  down  70^  of  the  aircraft  with  the  parameters  used. 

The  pajraimeters  used  are  not  specific  to  one  system,  but 
can  be  easily  changed  to  reflect  a specific  system. 

The  appendices  contain  results  of  the  programs  along 
with  the  calculations  used  to  determine  an  intercept  point 
for  the  aircraft  and  missile. 
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OPTICAL  THREAT  / LASER  COUNTERMEASURES 
ENGAGEMENT  ANALYSIS 

I . Introduction 

During  the  Arab-Israeli  War  in  1973f  the  use  of  opti- 
cally guided  surface-to-air  missiles  (SAMS)  and  anti-aircraft 
artillery  (AAA)  took  much  of  the  world  by  surprise.  It  was 
so  effective,  when  combined  with  conventional  systems,  that 
the  Arabs  were  able  to  seriously  challenge  Israeli  Air  supe- 
riority for  the  first  time. 

The  analysis  of  the  problems  faced  by  the  Israeli  Air 
Force  led  USAF  planners  to  the  conclusion  that  in  the  future,  ^ 

penetrating  aircraft  may  face  optically  and/or  electro-opti- 

Ij 

cally  directed  weapons.  Numerous  programs  have  since  been 
funded  to  develop  optical  receivers  and  countermeasures  de- 
vices. One  of  the  problems  encountered  was  the  need  for  a 
generic  computer  model  to  analyze  the  engagement  of  am  air- 
craft configured  with  the  proposed  optical  countermeasures 
equipment  against  threats  using  optically  directed  weapons. 

The  purpose  of  this  report  is  to  develop  such  a computer 
model. 

The  model  has  been  limited  to  am  engagement  of  one  air- 
craft against  one  threat.  Information  on  the  various  co\m- 
termeasures  devices  being  considered  or  developed  to  be  used 
in  this  role  was  not  available  when  this  report  was  written. 

Information  on  the  threats  against  which  these  devices  could 


1 


be  used  was  also  not  available  when  this  report  was  written. 
The  above  information  will  soon  be  available  at  the  Air  Force 
Avionics  Laboratory  Library  and  is  mostly  classified.  For 
these  reasons,  the  parameters  used  in  this  report  do  not  re- 
flect a specific  countermeasure  or  threat  device  but  rather 
a combination  of  generally  realistic  values.  Atmospheric 
propagation  has  been  ignored  due  to  time  constraints. 

The  initial  step  in  the  problem  was  an  extensive  liter- 
ature search.  Various  documents  were  obtained  to  aid  in  the 
preparation  of  this  report.  A Defense  Dociimentation  Control 
Center  search  was  initiated  to  obtain  technical  reports  per- 
taining to  the  problem.  Again,  these  documents  were  not 
available  or  their  classification  limited  their  use  to  gen- 
eralities and  not  specific  data. 

The  primary  problem  is  focused  around  the  computer 
model  itself.  It  was  required  that  this  model  be  flexible 
enough  to  allow  incorT)oration  at  a later  time  of  those  parts 
of  the  problem  not  covered  in  this  report  or  which  need  ex- 
panding. Also,  it  was  necessary  that  the  model  be  flexible 
enough  to  allow  changes  in  threat  and  OCM  parameters  to  sim- 
ulate any  of  the  threat  systems  and  proposed  receivers  and 
jamming  devices. 

For  these  reasons,  the  GASP  IV  simulation  language 
was  chosen  as  the  simulation  base.  It  combines  the  advan- 
tages of  being  written  in  FORTRAN,  immediately  available, 

Eind  modular  in  form.  Learning  to  use  GASP  IV  was  a major 
task  in  this  project. 
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II.  The  Scenario  and  GASP  IV 
Description  of  the  Scenario 

The  scenario  involves  one  aircraft  and  one  threat.  The 
purpose  is  to  analyze  their  engagement  by  collecting  statis- 
tics on  various  events  such  as  aircraft  entry  into  the  threat 
area  and  missile  firing.  The  final  goal  is  to  obtain  results 
on  aircraft  survivability.  All  distances  are  measured  in 
meters.  The  aircraft  initial  position  io  defined  by  POSAXI, 
POSAYI,  and  POSAZI  in  the  x,  y,  and  z directions,  with  the  z 
parameter  representing  aircraft  altitude.  The  threat  posi- 
tion is  defined  by  POSTX,  POSTY,  and  POSTZ.  Throughout  this 
simulation,  the  threat  has  been  located  at  the  origin. 

The  aircraft  flies  a straight  path  which  leads  directly 
over  the  threat.  Aircraft  maneuvers  were  not  incorporated 
due  to  time  constraints.  Also,  the  aircraft  velocity  and 
altitude  were  held  constant  in  each  simulation. 

The  threat  is  given  an  acquisition,  tracking,  and  fir- 
ing role.  Each  of  these  fimctions  has  limitations  which  are 
represented  by  a sphere  centered  at  the  threat.  The  outer 
surface  of  the  sphere  represents  the  limit  of  that  particular 
function.  As  an  example,  the  limit  of  an  acquisition  radar 
is  represented  in  Fig.  1 on  page  4,  The  aircraft  is  consid- 
ered to  be  in  range  of  the  acquisition  radar  when  it  enters 
this  sphere. 

Each  engagement  begins  with  the  aircraft  at  the  periph- 
ery of  acquisition.  The  aircraft  then  moves  through  the 
threat  area  and  passes  over  the  site.  When  the  aircraft 
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Fig.  1.  Effective  Range  of  Acquisition  Radar 


leaves  the  threat  sphere  of  influence  the  simulation  is 
ended. 

Fly  By 

In  planning  this  project,  it  was  decided  to  break  the 
problem  up  into  several  stages  to  be  accomplished  one  after 
the  other.  Each  of  these  pieces  is  a complete  program  in 
itself.  The  name  Fly  By  has  been  given  to  each,  followed  by 
a nximber  to  designate  where  it  appears  in  the  overall  pro- 
ject. Fly  By  1 is  the  initial  program.  Each  succeeding 
program  incorporates  the  preceding  ones  and  adds  a new  di- 
mension to  the  problem.  All  of  the  Fly  By  programs  are  de- 
scribed in  detail  in  Chapter  III. 
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GASP  IV 


To  speed  up  the  process  of  building  the  model,  a simu- 
lation language  was  needed  to  provide  structure  and  termi- 
nology. GASP  IV  is  such  a simulation  language.  The  programs 
built  using  GASP  IV  are  useful  as  both  models  and  analysis 
tools  (Ref  111). 

GASP  provides  the  user  with  a philosophy  and  a program- 
ming language.  Modeling  concepts  are  built-in  and  are  linked 
to  the  user  through  FORTRAN  routines  that  can  be  easily 
learned  and  used.  Other  advantages  are*  GASP  is  FORTRAN 
based  and  requires  no  separate  compiling  system?  it  is  mod- 
ular? it  is  easy  to  learn  because  it  is  FORTRAN  based?  it  is 
easily  modified  to  particular  applications?  and  it  can  be 
used  for  discrete,  continuous,  or  combined  simulation  (Ref 
1*3). 

GASP  uses  an  entity-attribute  approach  in  modeling.  En- 
tities are  those  things  in  a system  such  as  people,  equip- 
ment, and  raw  materials.  Attributes  are  the  characteristics 
of  these  entities.  Entities  having  a common  attribute  are 
grouped  into  files. 

In  simulating  a system,  one  desires  to  learn  something 
about  the  behavior  and  performance  characteristics  of  the 
system.  This  is  done  by  reproducing  the  activities  in  which 
the  entities  engage.  The  system  has  certain  states  defined 
in  terms  of  the  numeric  values  assigned  to  attributes  in  the 
system.  Simulation  causes  the  system  to  move  from  state  to 
state  and  is  the  dynamic  portrayal  of  the  states  of  a system 
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over  time.  The  two  most  important  activities  associated  with 
simulation  are  the  identification  of  the  entities  and  their 
attributes,  and  the  coding  of  the  attribute  values  into  char- 
acteristic system  states  (Ref  1j5)* 

A few  words  about  discrete,  continuous,  and  combined 
simulation  are  necessary.  As  a model  moves  from  state  to 
state,  entities  engage  in  activities  which  change  the  states. 
These  activities  are  marked  at  beginning  and  end  by  events. 
Take  as  an  example  an  elevator.  Consider  the  events  "eleva- 
tor starts  travel”  and  "elevator  arrives  at  next  floor”.  In 
a discrete  model,  the  time  of  arrival  of  the  elevator  at  the 
next  floor  would  be  known  from  the  event  "elevator  starts 
travel".  It  is  a given  value,  a discrete  value,  based  on 
the  time  at  which  the  elevator  started  moving.  In  a contin- 
uous model,  however,  the  time  between  floors  would  not  be 
deterministic.  The  functions  in  the  equation  of  motion  of 
the  elevator  would  have  variations  which  might  reflect  the 
nijmber  of  people  on  board,  the  direction  of  travel,  or  the 
number  of  serviceable  elevators.  The  time  between  floors 
would  be  variable  and  depend  on  these  and  other  factors. 

Such  attributes,  which  are  based  on  a continuous  representa- 
tion of  the  dynamics  of  the  attribute,  are  called  state  vari- 
ables in  GASP  (Ref  li8). 

When  the  dependent  variables  change  discreetly  at  spec- 
ified points  in  simulated  time,  discrete  simulation  is  oc- 
curring. When  they  change  continuously,  a continuous  simu- 
lation is  the  result.  And  when  the  dependent  variables 
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change  continuously  with  possible  discrete  jumps  superim- 
posed, the  result  is  defined  as  a combined  simulation  (Ref 
1.8). 

The  GASP  philosophy  is  explained  in  Chapter  2 of  refer- 
ence 1 . The  reader  will  also  find  printouts  and  flowcharts  j 

of  the  GASP  subroutines.  Numerous  examples  are  used  to  ex-  ! 

I 

plain  to  the  user  the  various  aspects  of  GASP  simulation. 

► 

[ The  most  important  of  the  GASP  subroutines  are  those 

f 

I written  by  the  user.  These  routines  will  be  mentioned  re- 

peatedly in  the  following  chapters.  Their  purposes  are  out- 
I lined  below. 

Subroutine  INTLC  is  used  to  initialize  non-GASP  vari- 
ables,  state  variables. 

I ables,  if  any.  SCOND  is  the  subroutine  which  determines  if 

a state-event  has  occurred.  When  either  a state  or  time- 
I event  occurs,  subroutine  EVNTS  is  used  to  determine  what 

[ action  to  take  for  each  event.  It  also  checks  to  see  which 

' state -event  has  occurred.  Subroutine  OTPUT  provides  a means 

for  the  user  to  obtain  output  in  addition  to  that  normally 
provided  by  GASP.  Two  other  subroutines  are  user  written 
but  were  not  used  in  this  project.  These  are  UERR,  which 
allows  specific  information  to  be  printed  for  errors,  and 

1 

[ SSAVE,  which  collects  data  for  the  user  to  print  out.  Dvanmy 

I 

I subroutines  were  used  for  UERR  and  SSAVE  throughout  the  sim- 

ulation. 

i 

GASP  Subroutines 

■ I 

I All  user  variables  have  been  listed  and  defined  on  | 

i ! 
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pages  v-viii.  Those  GASP  variables  necessary  for  under- 
standing the  discussions  in  this  report  are  listed  and  de- 
fined on  page  viii. 

There  are  three  GASP  subroutines  and  three  GASP  func- 
tions that  are  used  frequently  in  the  user  written  subrou- 
tines described  in  this  report.  They  are  defined  and  their 
functions  explained  in  the  following  paragraphs. 

Subroutine  COLCT(XX, ICLCT)  collects  sample  data  on  the 
variable  SS(XX).  It  calculates  the  mean,  standard  deviation 
of  the  mean,  standard  deviation,  coefficient  of  variation, 
minimiom,  maximum,  and  the  number  of  times  the  variable  was 
observed.  All  this  information  is  printed  out  in  tabular 
form  at  the  end  of  the  simulation  (Ref  1:133). 

Subroutine  HISTO(XX,IHIST)  tabulates  the  number  of 
times  a variable  is  within  a prescribed  cell  limit.  It  also 
calculates  observed,  relative,  and  cvimulative  frequency  for 
each  cell.  At  the  end  of  the  simulation,  it  prints  out  a 
plotted  histogram  of  the  relative  and  cumulative  frequency 
of  observations  (Ref  1:138). 

Each  of  these  subroutines  can  be  used  in  several  modes, 
depending  on  the  value  of  ICLCT  and  IHIST  used.  The  compu- 
tation and  reporting  mode  is  the  one  described  above  and  the 
only  one  used  in  this  report.  This  mode  is  indicated  by 
setting  positive  values  in  ICLCT  and  IHIST. 

The  other  subroutine  is  FILEM ( IFILE ) , which  is  used  to 
store  attributes  in  IFILE.  In  this  report,  FILEM  is  used  to 
change  the  attributes  of  time -events. 
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Two  of  the  functions  are  random  number  generators. 
RNORM(IPAR, ISTRM)  is  a normal  deviate  generator  which  uses 
the  random  number  stream  ISTRM  and  parameters  from  parameter 
set  IPAR.  BRAND (ISTRM)  uses  random  number  stream  ISTRM  to 
generate  a uniform  number. 

The  last  function  is  KROSS(IKRSG,IKRSD,CMULT,CADD,LDIR, 
TOL) . This  function  locates  and  specifies  state  conditions 
and  returns  a coded  value  indicating  whether  certain  condi- 
tions have  been  met.  IKRSG  is  the  index  of  the  crossing  var- 
iable. In  this  report,  it  will  always  be  a 1 , which  indi- 
cates that  the  crossing  variable  is  SS(1),  the  aircraft  range. 
There  is  no  crossed  variable,  so  IKRSD  will  always  be  0.  This 
value  is  multiplied  by  CMULT.  The  result  is  then  added  to 
CADD  to  obtain  a final  value  with  which  to  compare  the  cross- 
ing variable.  LDIR  and  TOL  indicate  the  direction  emd  tol- 
erance of  the  crossing  (Ref  1:121). 

When  f\mction  KROSS  is  called,  there  are  five  possible 
values  assigned  to  the  corresponding  LFIiAG.  A +2  value  in- 
dicates that  a positive  crossing  has  occurred  outside  the 
tolerance.  A +1  value  indicates  a positive  crossing  within 
tolerances,  while  a 0 indicates  no  crossing  at  all.  Values 
of  -2  and  -1  indicate  negative  crossings  outside  tolerance 
and  within  tolersince,  respectively. 
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III.  The  Fly  By  Programs 

Philosophy 

The  objectives  of  each  of  the  Fly  By  programs  is  dif- 
ferent. In  FBI,  the  primary  task  was  learning  to  use  GASP 

IV.  For  this  reason,  there  is  no  actual  detection  in  FBI. 

f 

The  aircraft  is  assumed  detected  as  soon  as  it  enters  the 
threat  range.  In  FB2,  FB3.  and  FB4,  however,  the  detection 
and  tracking  of  the  aircraft  becomes  probabilistic. 

The  primary  variable  in  determining  a probability  of 
detection  or  tracking  in  an  aircraft-ground  engagement  is 
the  aircraft  range  from  the  threat.  Thus,  SS(l)  is  used  in 
each  of  the  probabilities  computed.  But  there  are  other  pa- 
rameters which  affect  the  probabilities.  Some  of  these  are 
radar  power,  aircraft  signature,  and  weather  in  the  optics 
regime.  In  this  simulation,  these  parameters  are  not  in- 
cluded individually,  but  are  combined  into  one  number,  HIALT. 
Time  constraints  and  the  desire  to  keep  this  report  unclas- 
sified meant  that  HIALT  was  not  calculated,  but  simply  se- 
lected to  give  a realistic  probability  distribution  to  the 
detection  and  tracking  functions. 

In  determining  the  probability,  HIALT  is  divided  by 
range.  The  resulting  value,  PRBDET  or  PRBTRK,  is  then  com- 
pared with  1 . The  value  is  then  compared  with  a randomly 
selected  nvimber  from  a 0 to  1 uniform  distribution.  If  the 

1 

[ value  is  less  than  or  equal  to  the  raindom  number,  then  a yes 

i 

I condition  results.  This  means  the  aircraft  has  been  detect- 

I ed  or  that  the  tracking  device  has  locked  on. 


I 
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As  an  example,  assume  an  aircraft  at  TEFFR  of  a threat. 

Also  assume  the  following  values  applyi 

HI ALT= 5000.0 
TEFFR= 50000.0 

The  probability  of  detection  on  the  initial  scan  is  given 

HIALT/SS(1)=.10 

Thus,  there  is  a lOjS  probability  of  detection  on  the  first 

pass.  To  determine  if  detection  occurs,  a random  number  is 

selected.  If  it  falls  within  the  shaded  area  of  the  graph 

below,  then  the  aircraft  is  detected.  If  not,  detection 

does  not  occur,  and  the  program  continues,  I 


It  is  possible  for  the  aircraft  to  pass  the  site  and 
not  be  detected,  but  only  if  the  aircraft  altitude  is  great- 
er than  HIALT.  Assvime  the  aircraft  is  moving  500m/2s  scan, 
that  is  it  moves  500™  between  scans.  Then  the  probability 
of  it  not  being  detected  in  the  first  10  scans  with  the 
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above  variable  values  is  given  by 


10 

(l-HIALT/(SS(l)-500i)) 

i=0 

= 36.8;* 

In  FB2,  the  probability  of  detecting  an  aircraft  on  the 
first  scan  is  100^.  In  FB3  this  has  been  made  more  realis- 
tic by  changing  the  value  of  HIALT.  In  this  program,  the 
probability  of  detection  is  2.79%  on  the  first  scan.  This 
goes  to  100^  when  the  aircraft  is  5000m  from  the  threat.  In 
FB4,  the  probabilities  are  altered  so  that  detection  on  the 
first  scan  occurs  31.^  of  the  time.  The  probability  of 
lock-on  by  the  tracker  on  its  first  scan  is  95%> 

The  piirpose  of  these  distributions  is  not  to  show  actu- 
al detection  and  tracking  probabilities,  which  would  require 
the  radar  range  equation.  The  purpose  is  to  show  that  the 
program  is  capable  of  determining  the  probabilities  based 
on  the  data  input  of  the  user. 

Fly  By  1 

The  purpose  of  the  Fly  By  1 program  is  to  fly  an  air- 
craft through  a threat  area  on  a straight  path.  The  range 
from  the  aircraft  to  the  threat  is  computed  and  constantly 
updated.  When  the  aircraft  enters  the  effective  range  of 
the  acquisition  radar,  the  range  is  printed.  The  simulation 
ends  when  the  aircraft  passes  out  of  range  of  the  acquisition 
radar.  A general  flowchart  of  FBI  is  shown  in  Fig.  2 on 
page  13.  The  scenario  depicted  in  FBI  is  shown  in  Fig.  3 
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Fig.  2 A General  Flowchart  of  FBI 
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above . 

PROGRAM  FLYBYl  starts  the  simulation  by  setting  the 
cardreader  and  printer  numbers.  It  then  calls  GASP,  which 
initializes  GASP  variables  and  reads  in  the  data  cards.  From 
GASP,  the  program  goes  to  INTLC,  which  initializes  all  FBI 
variables  that  require  initialization.  A partial  listing  of 
INTLC,  showing  those  variables  initialized,  follows.  Only 
common  blocks  and  end  statements  will  be  omitted  from  the 
listings  found  in  this  report. 


P0SAXI=0.0 
POSAYI= 11 000.0 
P0SAZI=5000.0 
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POSAX=POSAXI 

POSAY=POSAYI 

POSAZ=POSAZI 

ASPEDX=0.0 

ASPEDY= -400.0 

ASPEDZ=0.0 

TEFFR= 10000.0 

P0STX=0.0 

P0STY=0.0 

P0STZ=0.0 


Once  initialization  has  been  accomplished,  INTLC  re- 
turns the  program  to  GASP,  which  then  calls  STATE.  In  STATE, 
TNOW  is  used  to  compute  the  aircraft  position  and  range. 

This  is  accomplished  with  the  following  statementst 


P0SAX=P0SAXI+ASPEDX*TN0W 
P0SAY=P0SAYI+ASPEDY*TN0W 
F0SAZ=P0SAZI+ASPEDZ*TN0W 
SS ( 1 ) = ( P0SAX**2+P0SAY**2+P0SAZ**2 ) ** . 5 


The  first  three  statements  take  the  initial  position  of  the 
aircraft  and  update  it  by  multiplying  the  speed  along  each 
direction  by  the  total  time  of  travel  in  that  direction.  The 
last  statement  then  computes  the  aircraft  range  from  the 
threat. 

On  return  from  STATE,  GASP  calls  SCOND,  which  checks  to 
see  if  a state  condition  has  occurred.  In  this  simulation, 
there  are  two  state-events.  These  are  aircraft  penetration 
into  the  threat  area  and  aircraft  withdrawal  from  the  threat 
area.  Checking  for  these  two  events  is  accomplished  in  SCOND 
by  the  following  statements i 


LFLAG(1)=KR0SS( 1,0, 0.0, 10000.0, +1,200.0) 
LFLAG(2)=KR0SS( 1,0, 0.0, 10000. 0,-1, 200.0) 
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In  the  first  statement,  the  aircraft  range  SS(1)  is  compared 
with  10000.0m  to  determine  if  a positive  crossing  within  a 
tolerance  of  200.0m  has  occurred.  When  this  crossing  occurs, 
LFLAG(l)  is  set  to  a value  of  +1,  indicating  that  the  air- 
craft is  leaving  the  threat  area.  In  the  second  statement, 
SS(1)  is  compared  with  10000.0m  to  see  if  a negative  cross- 
ing has  occurred  within  200.0m.  If  neither  state-event  has 
occurred,  GASP  updates  TNOW  and  returns  to  STATE  to  begin 
the  process  again.  This  continues  xmtil  one  of  the  state- 
events  occurs.  At  that  point,  GASP  calls  EVNTS. 

The  task  of  EVNTS  is  to  determine  which  state  or  time- 
event  has  occurred.  If  the  event  is  a time-event,  the  at- 
tribute ATRIB(2)  is  checked  to  see  which  of  the  time-events 
has  occurred.  There  are  no  time -events  in  FBI,  so  a de- 
tailed explanation  of  them  will  be  deferred  until  the  dis- 
cussion of  FB2.  For  state-events,  EVNTS  uses  the  following 
statements  to  determine  which  has  occurred  and  what  action 
to  take: 

IF  (LFLAG(1)=+1)  CALL  AORNG 

IF  (LFLAG(2)=-1)  CALL  AIRNG 

The  first  statement  checks  to  see  if  LFLAG(1)=+1.  If  yes, 
then  the  aircraft  is  out  of  range  and  AORNG  is  called.  The 
second  statement  checks  to  see  if  LFLAG(2)=-1,  indicating 
the  aircraft  is  within  range.  If  yes,  then  AIRNG  is  called. 
The  two  state-events  are  mutually  exclusive.  They  camnot 
occur  simultaneously  unless  an  error  has  been  made.  Also, 
when  EVNTS  is  called,  one  of  the  state-events  has  occurred 
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and  either  AORNG  or  AIRNG  must  be  called. 

Subroutine  AIRNG  is  called  when  the  aircraft  comes 
within  range  of  the  acquisition  radar.  It  computes  the  posi- 
tion of  the  aircraft  at  the  time  of  the  crossing  and  prints 
the  time,  coordinates,  and  range  of  the  crossing  event.  It 
does  this  with  the  following  statementsi 


POSAX=POSAXI+ASPEDX*TNOW 

POSAY=POSAyi+ASPEDy*TNOW 

POSAZ=POSAZI+ASPEDZ*TNOW 

PRIN  T ( 6 , 1 0 5 ) TNOW , POS AX , POSAY , POSAZ , SS ( 1 ) 

105  F0RMAT("TN0W  IS",F15.2  / "COORDINATES  ARE",3F15.2 
/ "RANGE  I 
1S",F15.2) 


Subroutine  AORNG  is  called  when  a positive  crossing 
occurs,  indicating  the  aircraft  has  passed  out  of  range  of 
the  radar.  The  time,  position,  and  range  of  this  event  are 
then  printed.  The  simulation  is  then  ended  by  resetting  the 
value  of  MSTOP  to  a negative  value.  The  statements  which 
accomplish  this  task  are  the  followingi 


PRINT (6,105)  TNOW , POSAX , POSAY , POSAZ , SS ( 1 ) 

105  FORMAT ("TNOW  IS",F15.2  / "COORDINATES  ARE",3F15.2 
/ "RANGE  I 
1S",F15.2) 

MST0P=-1 


Once  AORNG  has  been  called,  the  simulation  is  ended 
the  next  time  GASP  gains  control.  All  GASP  and  user  infor- 
mation requested  is  then  printed  out.  A sample  of  the  print- 
out from  FBI  is  listed  in  Appendix  A.  These  results  are  dis- 
cussed in  Chapter  V. 
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In  FBI,  the  aircraft  was  assumed  detected  when  it 
crossed  the  TEFFR,  In  FB2,  a more  realistic  detection  meth- 
od is  developed.  The  radar  is  allowed  to  begin  scanning  for 
the  aircraft  when  the  aircraft  crosses  the  TEFFR.  The  proba 
bility  of  the  aircraft  being  detected  by  the  radar  is  now  a 
function  of  range  to  the  aircraft  and  an  assigned  threat  ef- 
fectiveness number,  HIALT.  This  introduces  an  inverse  range 
randomness.  The  HIALT  variable  incorporates  all  the  param- 
eters which  would  normally  be  used  in  the  radar  range  equa- 
tion. For  FB2,  HIALT  was  chosen  to  give  a realistic  chance 
of  detection  at  TEFFR. 

Other  refinements  over  FBI  appear  in  the  printout. 
Headings  have  been  included  to  tell  when  the  aircraft  pene- 
trates, withdraws,  or  is  detected.  The  time,  coordinates, 
and  range  of  each  of  these  events  is  also  printed  out,  A 
general  flowchart  of  FB2  is  shown  in  Fig.  5 on  page  20.  In 
Fig.  4,  shown  on  page  19i  the  FB2  scenario  is  depicted. 

PROGRAM  FLYBy2  starts  the  simulation  by  setting  the 
cardreader  and  printer  nmbers.  It  then  calls  GASP,  which 
reads  in  the  data  cards  and  initializes  GASP  variables. 

GASP  then  calls  INTLC,  which  initializes  all  non-GASP  vari- 
ables. The  statements  which  accomplish  this  are  the  fol- 
lowings 

P0SAXI=0.0 
posAyi= 11 000.0 
P0SAZI=5000.0 
P0SAX=P0SAXI 
posAy=posAyi 


CALL  AIRNG 


POSAZ=POSAZI 

ASPEDX=0.0 

ASPEDY= -400.0 

ASPEDZ=0.0 

TANGLE = . 0349 

HIALT= 10000.0 

TEFFR=P0SAZ/STN ( TANGLE ) 

TEFFR=AMIN1 ( TEFFR , 1 0000 . 0) 

POSTX=0.0 

P0STY=0.0 

P0STZ=0.0 

SCANR=10.0 


As  can  be  seen,  there  are  a few  changes  over  FBI.  TANGI£  is 


% 


Fig.  5.  A General  Flowchart  of  FB2. 
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the  terrain  masking  angle  in  radians.  It  is  used  to  deter- 
mine how  far  above  the  horizon  the  radar  can  "see",  assuming 
a flat  earth.  The  threat  effective  range  is  calculated  from 
this  angle  and  then  set  equal  to  the  minimum  of  the  actual 
range  of  the  radar  and  its  range  at  the  aircraft  altitude  to 
be  considered  in  this  simulation.  Another  change  is  SCANR, 
which  represents  the  sweep  rate  of  the  radar  in  seconds. 

Once  initialization  has  been  accomplished,  INTLC  re- 
turns to  GASP,  which  then  calls  STATE.  In  STATE,  as  in  FBI, 
TNOW  is  used  to  compute  the  aircraft  position  and  range. 

The  statements  used  to  accomplish  this  are  the  same  used  in 
FBI.  They  are  listed  on  page  15. 

On  return  from  STATE,  GASP  calls  SCOND  and  a check  is 
made  for  state  conditions.  As  in  FBI,  aircraft  penetration 
and  aircraft  withdrawal  are  the  only  state-events.  The  fol- 
lowing statements  are  used  in  SCOND  to  check  for  these 
events! 

LFLAG(1)=KR0SS(1,0,0.0,TEFFR,+1, 200.0) 

LFLAG(2)=KROSS(1,0,0.0,TEFFR, -1,200.0) 

In  the  first  statement,  SS(1)  is  compared  to  the  variable 
TEFFR  to  determine  if  a positive  crossing  has  occurred.  This 
would  indicate  aircraft  withdrawal.  In  statement  two,  a 
check  for  aircraft  penetration  is  made  by  comparing  SS(1) 
and  TEFFR  for  a negative  crossing.  The  tolerance  in  both 
cases  is  200.0m,  the  same  as  in  FBI.  If  neither  state-event 
has  occurred,  GASP  updates  TNOW  and  returns  to  STATE.  This 
process  continues  \jntil  an  event  is  detected  in  SCOND  and 
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GASP  calls  EVNTS. 

EVNTS  determines  which  state  or  time -event  has  occurred. 
The  argiJiment  IX  in  subroutine  EVNTS (IX)  is  used  to  distin- 
guish time  and  state-events.  All  time-events  are  initialized 
with  at  least  two  attributes.  ATRIB(l)  is  the  time  of  occur- 
rence of  the  event,  while  ATRIB(2)  is  the  event-code  of  that 
time -event.  Initialization  also  includes  the  setting  of  the 
state-event  code,  IIEVT,  which  is  unique  from  all  ATRIB(2) 
time-event  identifiers.  The  IX  argument  is  set  by  GASP.  If 
TNOW  corresponds  to  a time-event  occurrence  time,  ATRIB(l), 
then  IX  is  set  to  the  code  for  that  event,  ATRIB(2).  If 
SCOND  detects  a state-event,  GASP  sets  IX  to  IIEVT.  A par- 
tial listing  of  EVNTS  appears  belowt 

GO  TO  (101, 102), IX 
102  CALL  SCAN 
GO  TO  110 
101  CONTINUE 

IF  (LFLAG(2) .EQ.-l)  CALL  AIRNG 
IF  (LFLAG(l) .EQ.+l)  CALL  AORNG 
110  RETURN 

The  computed  GO  TO  statement  separates  time -events  from 
state-events.  In  FB2,  the  state-event  code  is  a 1 while  the 
only  time-event  is  identified  by  2.  The  time  of  occurrence 
of  a time-event  is  determined  by  ATRIB(l)  of  the  time-event. 
This  attribute  is  initially  set  on  data  cards  and  can  be 
changed  during  the  simulation.  When  TN0W=ATRIB(1 ) of  any 
time-event,  GASP  sets  IX-ATRIB(2)  of  the  same  time-event. 

When  a state-event  is  detected  in  SCOND,  GASP  sets  IX  to  the 
state-event  code. 


22 


If  the  event  is  a time-event,  IX=2  causes  EVNTS  to  call 


SCAN.  In  this  case,  there  is  no  need  to  check  for  state- 
event  flags,  so  a return  to  GASP  is  initiated  after  return 
from  SCAN.  If  the  event  is  a state -event,  then  all  flags 
are  checked  to  see  which  event  has  occurred  and  what  action 
to  take.  For  aircraft  penetration,  LFLAG(2)=-1  and  AIRNG  is 
called.  For  aircraft  withdrawal,  LFIiAG(l)=+l  and  AORNG  is 
called.  The  three  subroutines  which  EVNTS  can  call  are  dis- 
cussed in  the  following  order:  AORNG,  AIRNG,  and  SCAN. 

As  noted,  AORNG  is  called  for  aircraft  withdrawal.  The 
time,  position,  and  range  of  aircraft  withdrawal  are  printed 
and  the  simulation  ended  using  the  statements  listed  below. 


PRINT*,"  THE  AIRCRAFT  IS  OUT  OF  RANGE  " 

PRINT (6,105)  TNOW , POSAX , POSAY , POSAZ , SS ( 1 ) 

105  FORMAT ("  TNOW  IS  ",F15.2  / " COORDINATES  ARE  ", 
3F15.2  / " RANGE  I 
IS  ",F15.2) 

MST0P=-1 


Subroutine  AIRNG  has  a few  changes  from  FBI  due  to  the 
time -event.  A partial  listing  of  AIRNG  appears  belowt 


PRINT*,"  THE  AIRCRAFT  IS  IN  RANGE  " 

PRINT ( 6,105)  TNOW , POSAX , POSAY , POSAZ , SS ( 1 ) 

105  FORMAT ("  TNOW  IS  ",F15.2  / " COORDINATES  ARE  ", 
3FI5.2  / " RANGE  I 
IS  ",F15.2) 

ATRIB(1)=TN0W 
ATRIB(2)=2.0 
CALL  FILEM(l) 


The  first  print  statement  just  tells  the  user  that  the  air- 
craft is  in  range.  The  time,  position,  and  range  of  the 

event  are  then  printed.  In  addition,  since  the  aircraft  is 

i 

I 
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in  range,  the  radar  can  start  scanning.  Thus,  the  time  of 
the  event  SCAN  is  set  to  TNOW  and  its  identifier  to  2.0. 

GASP  subroutine  FILEIVl(l)  is  used  to  file  the  changed  attri- 
butes in  file  1.  Immediately  upon  return  from  AIRNG,  GASP 
notes  that  ATRIB(l)  for  time -event  number  2 is  equal  to  TNOW. 
It  then  calls  EVNTS  with  IX=2  and  EVNTS  calls  SCAN. 

The  purpose  of  SCAN  is  to  look  at  the  aircraft  at  spec- 
ified intervals  and  determine  if  the  aircraft  has  been  de- 
tected. The  statements  used  to  accomplish  this  are  the  fol- 
lowing! 

PRBDET=AIVIIN1  (HIALT/SS  ( 1 ) , 1 . ) 

X=DRAND(1) 

IF  (X.LE.PRBDET)  CALL  DETECT 
ATRIBd  ) = TNOW+SCANR 
ATRIB(2)=2.0 
CALL  FILEM(l) 

The  first  statement  determines  the  probability  of  detection 
by  dividing  HIALT  by  SS(1)  and  taking  the  smallest  of  this 
value  and  1.  The  purpose  of  the  next  statement  is  to  take  a 
random  nvimber  X from  the  GASP  generator  BRAND.  X is  then 
compared  with  the  PRBDET  number.  If  it  is  greater  than 
PRBDET,  then  detection  has  not  occurred.  The  time  of  event 
SCAN,  ATRIB(l),  is  set  to  TNOW+SCANR,  which  means  the  radar 
must  sweep  before  trying  to  detect  again.  ATRIB(2)  is  also 
reset  and  both  values  put  in  file  1 by  FILEM. 

If  X is  less  than  or  equal  to  PRBDET,  then  detection 
has  occurred  and  DETECT  is  called.  A listing  of  the  state- 
ments used  in  DETECT  follows i 
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■ 1 

I 

PRINT*,"  THE  AIRCRAFT  HAS  BEEN  DETECTED  " i 

PRINT (6, 105)  TNOW.POSAX.POSAY.POSAZ.SSd)  I 

105  FORMAT ("  TNOW  IS  ",F15.2  / " COORDINATES  ARE  ",  ! 

3F15.2  / " RANGE  I 1 

IS  ",F15.2) 

MSTOP=-l 

The  fact  that  the  aircraft  has  been  detected  is  first  printed, 
and  then  the  time,  position,  and  range  of  the  detection.  The 
simulation  is  ended  by  setting  MST0P=-1. 

It  is  important  to  note  that  if  the  aircraft  is  detect- 
ed, it  is  not  allowed  to  continue  out  of  range.  Either 
DETECT  or  AORNG  can  be  called  in  FB2,  but  never  both.  Re- 
sults of  the  FB2  program  are  listed  in  Appendix  A and  dis- 
cussed in  Chapter  V. 

Fly  By  2 

The  purpose  of  the  Fly  By  3 program  is  to  make  a speci- 
fied niamber  of  rnns  with  FB2,  collect  statistics  on  each  run, 
and  then  print  the  final  results  of  all  the  r\ms.  It  was 
determined  that  the  way  to  do  this  was  to  reinitialize  each 
run  at  the  point  where  FB2  had  been  ended.  The  covinter 
NACFT  is  set  in  PROGRAM  FLYBY3  to  determine  how  many  runs 
will  be  made.  Each  loin  is  a FB2  program,  with  certain 
changes  to  reflect  more  realistic  parameters,  collect  sta- 
tistics, and  to  continue  or  stop  according  to  the  counter 
NACFT.  For  this  reason,  only  those  changes  made  to  FB2  will 
be  noted  2ind  explained.  All  other  subroutines  and  statements 
are  the  same  as  in  FB2.  A general  flowchart  of  FB3  is  shown 
in  Fig.  6 on  page  26. 

There  are  few  chainges  to  INTLC.  Most  of  these  reflect 
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changes  to  more  realistic  numbers.  These  new  statements 
are  shewn  belowi 


TN0W=0.0 

HIALT=4000.0 

P0SAYI=42000.0 

ASPEDY=- 200.0 

TEFFR= AMINl ( TEFFR , 40000 . 0 ) 

SCANR=2.0 


The  first  statement  is  crucial  in  that  it  resets  the  time 
when  INTLC  resets  the  aircraft  to  its  initial  position. 

AORNG  is  the  next  subroutine  which  has  been  changed 
from  FB2.  A listing  of  AORNG  is  shown  below* 


PRINT*,"  THE  AIRCRAFT  IS  OUT  OF  RANGE 

ATRIB(1)=1.E20 

ATRIB(2)=2.0 

CALL  FILEM(l) 

NACFT=NACFT-1 
IF  (NACFT.EQ.O)  GO  TO  110 
CALL  INTLC 
GO  TO  111 

110  MST0P=-1 

111  RETURN 


First,  the  flag  " THE  AIRCRAFT  IS  OUT  OF  RANGE  " is  printed. 
This  means  that  the  aircraft  cbxi  no  longer  be  detected.  To 
stop  the  detection  process,  the  time  of  occurrence  of  SCAN 
is  set  far  into  the  future.  Next,  the  number  of  aircraft  is 
reduced  by  1 . If  there  are  no  aircraft  remaining,  the  simu- 
lation is  ended  by  setting  MST0P=-1.  If  the  counter  NACFT 
is  not  0,  however,  AORNG  calls  INTLC  to  begin  the  simulation 
over  with  a new  aircraft.  No  statistics  are  collected  on 
aircraft  withdrawals  in  FB3. 

The  only  other  subroutine  which  is  different  from  FB2 
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is  DETECT.  A listing  of  this  subroutine  follows i 


PRINT*,"  THE  AIRCRAFT  HAS  BEEN  DETECTED  " 

PRINT( 6 ,105)  TNOW , POSAX , POSAY , POSAZ , SS ( 1 ) 

105  FORMAT("  TNOW  IS  ",F15.2  / " COORDINATES  ARE  ", 
3FI5.2  / " RANGE  I 
IS  ",F15.2) 

ATRIB(1)=1.E20 
ATRIB(2)=2.0 
CALL  FILEM(l) 

CALL  HISTO  (SS(1),1) 

CALL  COLCT  (SS(1),1) 

NACFT=NACFT-1 
IF  (NACFT.EQ.O)  GO  TO  700 
CALL  INTLC 
GO  TO  110 
700  MST0P=-1 
110  RETURN 


The  fact  that  the  aircraft  has  been  detected  is  first  noted 
by  a print  statement,  and  then  the  time,  position,  and  ramge 
of  detection  is  printed.  Since  detection  is  as  far  as  the 
program  goes  with  each  aircraft,  the  next  step  is  to  stop 
scanning  this  aircraft,  reinitialize,  and  go  on  to  another 
aircraft,  if  any.  To  stop  the  SCAN  subroutine  from  being 
called  again,  its  time  of  next  occurrence  is  set  far  into 
the  futvire . HISTO  aind  COLCT  are  then  used  to  collect  sta-  • 
tistics  on  the  range  of  detection.  The  number  of  aircraift 
is  reduced  by  1 and  if  any  remain,  a call  to  INTLC  starts 
the  simulation  over.  If  there  are  no  more  aircraft,  MSTOP 
is  set  to  a negative  value  to  stop  the  simulation. 

When  the  simulation  is  ended,  GASP  prints  out  the  re- 
quired data.  In  FB3f  COLCT  will  give  various  parameters  of 
the  range  of  detection,  such  as  minimimi  aind  maximum  values. 
HISTO  prints  out  a histogram  showing  where  the  aircraft  was 
detected  on  each  mm.  A discussion  of  the  results  of  FB3 
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can  be  found  in  Chapter  V.  A sample  of  these  results  is 
given  in  Appendix  A. 


Fly  By  4 

Fly  By  4 is  the  last,  aind  therefore  the  most  complicated, 
of  the  Fly  By  programs.  There  are  several  additions  over  its 
predecessor,  FB3.  Most  importamt  of  these  is  the  addition 
of  another  radar  to  the  threat.  This  is  a tracking  device. 

It  tracks  the  aircraft,  computes  its  position,  and  then  calls 
subroutine  FIREM  to  fire  missiles  at  the  target.  In  addition 
to  collecting  statistics  on  range  of  detection,  the  program 
also  collects  statistics  on  the  rauige  at  which  the  tracking 
device  "locks -on"  and  the  rainge  at  which  a missile  is  fired 
at  the  aircraft.  Other  output  includes  the  number  of  sur- 
viving aircraft  and  the  number  shot  down.  The  FB4  scenario 
is  shown  in  Fig.  7 on  page  30.  A general  flowchart  of  FB4 
is  shown  in  Fig.  8 on  pages  32-34. 

PROGRAM  FLYBY4  starts  off  the  simulation  with  the  fol- 
lowing statements I 

NCRDR=5 
NPRNT=6 
NACFT=  20 
NACFTI=NACFT 
NSURR=NACFT 
CALL  GASP 
CALL  EXIT 

The  cardreader  and  printer  numbers  are  set  first,  and  then 
the  counters  NACFT,  NACFTI,  and  NSURR.  The  main  program  then 
calls  GASP. 

After  initializing  GASP  variables  and  reading  in  all 
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CALL  AIRNG 
(initiate  scan) 


CONTINUE 

SCANV 


- CALL  DETECT 
(when  acft  is 
detected) 


TEFFR 


CALL  AIRNGT 
(initiate 
tracker  scan) 


CONTINUE 

SCAN-y 


CALL 

TRACK  (when 
lock-on  occurs) 


EFFRT 


CALL  FIREM  - ^ ^ 


N 

EFFRM 


_ - CALL  AORNG 
^ ( end ) 


MINERF - ^ 


(0,0,0) 


Fig.  7.  FB4  scenario 


data  cards,  GASP  calls  INTLC.  A partial  listing  of  the 
statements  in  this  subroutine  follows i 


MFLAG=1 

KFLAG=0 

JFLAG=1 


EFFRT= 50000.0 

MINEFR= 14000.0 

EFFRM=40000.0 

TN0W=0.0 

POSAXI=0.0 

P0SAYI=401000.0 

P0SAZI=5000.0 

P0SAX=P0SAXI 

P0SAY=P0SAYI 

P0SAZ=P0SAZI 

ASPEDX=0.0 

ASPEDY=- 222.0 

ASPEDZ=0.0 

HIALT=4000.0 

TANGLE= . 0349 

TEFFR=POSAZ/SIN ( TANGLE ) 

TEFFR= AMIN 1(TEFFR, 400000.0) 

P0STX=0.0 

P0STY=0.0 

P0STZ=0.0 

SCANR=2.0 

SCANRT= . 5 

NMISL=2 

DELAY=0.0 

RVEL=-222.0 

XMSPED=444.0 


The  three  flags  MFLAG,  KFLAG,  and  JFLAG  are  used  to  deter- 
mine if  certain  conditions  have  been  met.  They  will  be  de- 
scribed in  detail  in  the  descriptions  of  the  subroutines 
where  they  are  used.  The  effective  range  of  the  tracking 
device,  minimum  effective  range  of  the  missile,  and  the  mis- 
siles maximum  effective  range  are  read  in  as  EFFRT,  MINEFR, 
and  EFFRM,  respectively.  Otrer  new  parameters  are  DELAY,  a 
timing  device  used  in  FIREMi  NMISL,  the  number  of  missiles 
the  threat  is  allowed  to  fire  at  each  aircraft;  RVEL,  the 
speed  of  the  aircraft;  and  XMSPED,  the  missile  speed.  After 
initialization,  INTLC  returns  to  GASP,  which  then  calls 
STATE. 

Subroutine  STATE  is  used  in  FB4  exactly  the  same  as  in 
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Fig.  8.  A General  Flowchart  of  FB4  (Concluded) 
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FBI,  FB2,  and  FB3.  It  updates  the  aircraft  position  and  the 
aircraft  range  SS(1),  the  only  state  variable.  STATE  is 
listed  on  page  15. 


Upon  exiting  from  STATE,  GASP  calls  SGOND,  which  deter- 
mines if  any  state  conditions  have  been  met.  The  statements 
in  SCOND  are  as  follows  1 

LFLAG(1)=KR0SS(1,0,0.0,TEFFR, +1,200.0) 

LFLAG(2)=KR0SS(1,0,0.0,TEFFR, -1,200.0) 

lflag(3)=kross(i,o,o.o,effrt, -1,50.0) 

LFLAG (4)=KR0SS (1,0, 0.0, MINE FR, -1,50.0) 

LFLAG(l)  and  LFLAG (2)  are  used  the  same  way  as  in  previous 
programs.  They  check  to  see  if  the  aircraft  is  withdrawing 
or  penetrating,  respectively.  LFLAG(3)  and  LFLAG(4)  reflect 
the  addition  of  two  new  state  conditions.  When  the  aircraft 
comes  within  effective  range  of  the  tracking  device,  LFLAG(3) 
will  be  set  to  a value  of  -1  to  denote  the  crossing.  When 
LFLAG(4)  is  set  to  -1,  the  aircraft  has  passed  witnin  mini- 
mum effective  range  of  the  missile  and  cannot  be  fired  on. 

It  is  assumed  in  this  scenario  that  missiles  will  be  fired 
at  incoming  aircraft  only. 

The  program  progresses  with  GASP  constantly  updating 
TNOW  until  one  of  the  events  in  SCOND  occurs.  A call  to 
EVNTS  is  then  made  with  the  appropriate  IX  value.  The  first 
call  to  EVNTS  must  be  a state-event  since  both  time -events 
have  ATRIB(l)  read  in  as  1.E20.  The  statements  in  EVNTS 
are  as  follows i 

GO  TO  (101, 102, 103), IX 
102  CALL  SCAN 
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GO  TO  110 
103  CALL  SCANT 
GO  TO  110 
101  CONTINUE 

IF  (LFLAG(2) .EQ.-l)  CALL  AIRNG 
IF  (LFLAG(l) .EQ.+l)  CALL  AORNG 
IF  (LFLAG(3) .EQ.-l)  CALL  AIRNGT 
IF  (LFLAG(4) .EQ.-l)  CALL  OTPUT 

The  computed  GO  TO  statement  separates  time-events  from 
state-events.  In  FB4,  the  state-event  code  is  1.  If  IX=1, 
then  each  flag  is  checked  and  the  appropriate  subroutine  is 
called.  Each  of  the  subroutines  will  be  discussed  in  the 
following  paragraphs. 

The  AIRNG  subroutine  is  the  same  as  in  FB3  except  for 
the  print  statements.  A listing  of  this  subroutine  followst 

ATRIB(1)=TN0W 
ATRIB(2)=2.0 
CALL  FILEM(l) 

The  purpose  of  this  subroutine  is  to  initiate  the  scanning 
of  the  acquisition  radar.  This  is  a time-event,  SCAN,  and 
it  is  started  by  setting  the  start  time  to  TNOW.  ATRIB(2) 
is  used  to  identify  which  time-event  is  being  changed. 

FILEM  does  the  storage  and  retrieval  job. 

After  return  to  EVNTS  from  AIRNG,  EVNTS  gives  the  pro- 
gram back  to  GASP.  GASP  notes  that  ATRIB(l)  of  time- event  2 
is  equal  to  TNOW.  SCAN  is  called  immediately.  A listing  of 
this  subroutine  follows i 


PRBDET=AMIN1 (HIALT/SS( 1 ) , 1 . ) 
X=DRAND ( 1 ) 

IF  (X-PRBDET)101,101,102 
101  CALL  DETECT 
GO  TO  103 


r ■ 1 


102  KFLAG=1 

103  ATRIB(1)=TN0W+SCANR 
ATRIB(2)=2.0 

CALL  FILEM(l) 

110  RETURN 


The  first  statement  determines  the  probability  of  detecting 
the  aircraft,  PRBDET,  Next,  a random  number  X is  called 
from  BRAND  and  compared  with  PRBDET,  If  X is  greater  than 
PRBDET,  detection  does  not  occur.  The  flag  KFLAG  is  then 
set  to  1.  This  means  if  the  aircraft  is  detected  and  then 
lost,  the  fact  that  detection  has  not  been  continuous  will 
be  stored  in  KFLAG.  The  last  statements  reassign  a new  time 
for  SCAN  to  begin  again,  depending  on  the  scan  rate  of  the 
radar , SCANR . 

If  X is  less  than  PRBDET,  then  the  aircraft  has  been 
detected  and  DETECT  is  called.  A listing  of  this  subrou- 
tine appears  be low i 


IF  (MFLAG.EQ.O)  GO  TO  800 
CALL  HISTO  (SS(1),1) 

CALL  COLCT  (SS(1),1) 

MFLAG=0 

TDETI=TNOW 

800  IF  (KFLAG .EQ.l)  GO  TO  900 
TTIME=TNOW-TDETI 
GO  TO  110 
900  TDETI=TNOW 
KFLAG =0 
110  RETURN 


The  first  statement  uses  MFLAG  to  determine  if  this  is  the 
first  time  the  aircraft  has  been  detected.  If  yes,  then 
statistics  are  collected  on  range  of  detection.  MFLAG  is 
then  changed  so  that  statistics  on  range  of  detection  can- 
not be  collected  on  this  aircraft  again.  Also,  the  initial 
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time  of  detection,  TDETI,  is  set. 


{ 


I 


i 
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If  detection  was  not  made  during  the  previous  scan,  the 
initial  time  of  detection  is  reset.  KFLAG  is  used  to  deter- 
mine if  detection  did  not  occur  and  TDETI  is  reset  only  if 
KFLAG=1.  If  TDETI  is  reset,  KFLAG  is  also  reset.  If  detec- 
tion has  been  continuous,  the  total  time  is  computed  and 
stored  in  TTIME.  This  process  of  scanning  and  detecting 
continues  until  another  state-event  or  time-event  occurs. 

In  FB4,  this  happens  when  the  aircraft  crosses  EFFRT  and 
comes  within  range  of  the  tracking  device.  AIRNGT  is  then 
called  by  EVNTS. 

The  AIRNGT  subroutine  serves  the  same  purpose  for  the 
tracking  device  as  AIRNG  does  for  the  acquisition  radar.  A 
listing  of  this  subroutine  appears  belowi 

ATRIB(1)=TN0W 
ATRIB(2)=3.0 
CALL  FILEM(l) 

The  time  of  occurrence  of  time -event  3i  which  is  SCANT,  is 
set  to  TNOW.  Upon  return  to  GASP,  SCANT  is  called  immedi- 
ately. 

SCANT  is  to  the  tracking  device  what  SCAN  is  to  the 
acquisition  radar.  It  scans  the  aircraft  at  periodic  in- 
tervals and  determines  if  the  aircraft  is  being  tracked  by 
the  tracking  device.  The  following  is  a listing  of  SCANTi 


IF  (TTIME-10.0)  100,100,101 

100  PRBTRK=0.0 
GO  TO  102 

101  PRBTRK=AMIN1(HIALT/SS(1),1.) 
Y=RN0RM(1,1) 
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102  IF  (y-PRBTRK)  103,103,104 

103  IF  (JFLAG.EQ.O)  GO  TO  105 
JFLAG=0 

CALL  COLCT  (SS(l),2) 

CALL  HISTO  (SS(1),2) 

105  CALL  TRACK 

104  ATRIB(1)=TN0W+SCANRT 
ATRIB(2)=3.0 

CALL  FILEM(l) 

110  RETURN 


The  first  statement  determines  if  the  aircraft  has  been  de- 
tected long  enough  for  the  tracking  device  to  get  the  infor- 
mation necessary  for  accurate  positioning  to  find  the  air- 
craft. An  arbitrary  time  of  10.0s  has  been  used.  If  TTIME 
is  not  greater  than  10.0s,  "lock-on"  is  not  allowed.  If 
TTIME  is  greater  than  10.0s,  then  the  probability  of  track- 
ing, PRBTRK,  is  computed  in  the  same  manner  as  PRBDET  in 
SCAN.  A random  number  Y is  then  pulled  from  BRAND  and  if 
it  is  greater  than  PRBTRK,  the  aircraft  is  not  detected  by 
the  tracking  device.  When  this  occurs,  the  time  of  next 
occurrence  of  SCANT  is  set  to  TNOW+SCANRT,  where  SCANRT  is 
the  scanning  rate  of  the  tracking  device. 

If  Y is  less  than  or  equal  to  PRBTRK,  the  tracking  de- 
vice ”locks-on"  to  the  aircraft.  The  flag  JFLAG  is  then 
used  to  determine  if  this  is  initial  "lock-on".  If  yes, 
then  statistics  are  collected  and  JFLAG  set  so  that  they 
cannot  be  collected  again  on  this  aircraft.  Then  subrou- 
tine TRACK  is  called.  If  statistics  have  already  been  col- 
lected, TRACK  is  called  immediately. 

TRACK  is  used  to  determine  where  the  aircraft  is  and 
to  compute  an  intercept  point  at  which  to  fire  the  missile. 
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A listing  of  TRACK  follows i 


D=RN0RM(2,2) 

E=RN0RM(3,3) 

F=RN0RM(4,4) 

SST=SS(1)+D 

THETA=ATAN(POSAZ/POSAY)+E 

PHE=ATAN(P0SAY/P0SAX)+F  ' 

BARK= ( (POSAY-POSTY )/(P0SAZ-P0STZ ) )**2. 0 

BITE=BARK*XMSPED**2. 0 

D0G1=BITE+XIVISPED**2.  0 

D0G2=-BITE+RVEL**2.0 

D0G3=4.0*(RVEL*XMSPED)**2.0 

D0G4=D0G3**.5 

ANGLE=ACOS( ( -D0G4+ (D0G3-4*D0G1*D0G2)**. 5)/( 2*D0G1 ) ) 

MFLTT= (POSAZ-POSTZ )/{ -RVEL*SIN (ANGLE ) ) 

POSAX=POSAXI+ASPEDX* ( TNOW+MFLTT) 

POSAY=POSAYI+ASPEDY*(TNOW+MFLTT) 

POSAZ=POSAZI+ASPEDZ*( TNOW+MFLTT) 

SST=(P0SAX**2+P0SAY**2+P0SAZ**2)**.5+D 

PHE=ATAN (POSAY/POSAX )+F 

THE  TA= ATAN ( POSAZ/POSAY )+E 

POSMX=SST*COS(PHE) 

POSMY=SST*COS (THETA) 

POSMZ=SST*SIN (THETA ) 

IF  (SS(1)-EFFRM)101,101,110 
101  CALL  FIREM 
110  RETURN 


D,  E,  and  F are  random  numbers  pulled  from  RNORM.  They  are 
then  added  to  SST,  THETA,  and  PHE,  respectively,  to  simulate 
errors  in  the  tracking  device.  The  range,  elevation,  and 
azimuth  of  the  aircraft  as  seen  by  the  tracker  are  then  used 
along  with  aircraft  and  missile  speed  to  determine  an  inter- 
cept point.  The  formulas  and  diagrams  used  to  compute  this 
intercept  are  listed  in  Appendix  B.  The  quadratic  equation 

-b  ± (b^  - 4ac)^/^ 


is  used  to  make  the  final  calculation.  The  BARK,  BITE,  and 
DOG  statements  are  all  various  pieces  of  this  equation.  The 
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final  answer  is  the  time  of  flight  of  the  missile  to  the 
intercept  point,  or  MFLTT. 


Once  MFLTT  is  known,  the  computed  aircraft  position  is 
updated  by  the  appropriate  time.  A new  range,  elevation, 
and  azimuth  are  then  computed  and  used  to  determine  a posi- 
tion at  which  to  fire  the  missile.  These  coordinates  are 
POSMX,  POSMY,  and  POSMZ . The  missile,  when  fired,  is  as- 
sumed to  make  a perfect  launch  and  flight  to  these  coordi- 
nates and  detonate. 

To  determine  if  a missile  is  to  be  fired,  the  range  to 
the  aircraft  is  compared  with  the  missile's  effective  range. 
If  the  aircraft  is  not  in  range,  the  program  is  returned  to 
SCANT  where  the  process  continues.  When  the  aircraft  is  in 
range,  TRACK  calls  FIREM. 

The  FIREM  subroutine  is  used  to  fire  missiles  at  the 
aircraft.  A listing  of  FIREM  follows i 


IF  (NMISL.GT.O)  GO  TO  108 
GO  TO  111 

108  IF  (TNOW.GE. DELAY)  GO  TO  109 
GO  TO  110 

109  XMISS=((P0SMY-P0SAY)**2.0+ (POSMX -POSAX)"*''  0+ 

(POSMZ -POSAZ ) **2 . 0 ) ** 

1.5 

DELAY=TN0W+5.0 

NMISL=NMISL-1 

CALL  C0LCT(SS(1),3) 

CALL  HIST0(SS(1),3) 

IF  (XMISS-66. 0)101, 110, 110 
101  NSURR=NSURR-1 
111  CALL  OTPUT 

110  RETURN 


A check  is  first  made  to  determine  if  the  missiles  at  the 
site  have  all  been  fired.  If  yes,  then  an  immediate  call 
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is  made  to  OTPUT.  If  no,  then  the  time  since  the  last  mis- 
sile was  fired  is  computed  to  determine  if  another  missile 
is  ready.  The  time  DELAY  is  used  to  prepare  a missile  for 
firing  after  a missile  is  fired.  If  another  missile  is  not 
ready,  a return  to  TRACK  is  made. 

If  a missile  is  ready  for  firing,  it  is  fired  by  com- 
puting the  miss  distance  XMISS.  Statistics  are  then  col- 
lected on  the  range  at  which  firing  was  accomplished.  Also, 
the  number  of  missiles  is  reduced  by  1 and  the  time  delay 
between  missiles  is  set.  To  determine  if  the  aircraft  has 
been  shot  down,  the  distance  missed  is  compared  with  66m, 

The  aircraft  is  considered  shot  down  if  XMISS  is  less  than 
66m.  If  the  aircraft  is  not  shot  down,  the  program  returns 
to  TRACK.  If  the  aircraft  is  shot  down,  then  the  number  of 
survivors  is  reduced  by  1 and  OTPUT  is  called. 

A listing  of  subroutine  OTPUT  appears  belowi 

NNSURR=NACFTI -NSURR 
NACFT=NACFT-1 
IF  (NACFT.EQ.O)  GO  TO  200 
CALL  INTLC 
GO  TO  110 

200  PRINT  (6,1 60)  NSURR, NNSURR 

160  FORMAT  ("  THE  NUMBER  OF  SURVIVING  ACFT  IS  ",I5./i 
" THE  NUMBER  SHOT 
1 DOWN  IS  ”,15) 

700  MST0P=-1 
110  RETURN 

This  subroutine  first  computes  the  nxmber  of  aircraft  shot 
down,  NNSURR.  The  number  of  aircraft  is  reduced  by  1 and 
then  OTPUT  determines  if  any  aircraft  remain  to  fly  through 
the  threat  area.  If  yes,  INTLC  is  called  and  the  simulation 
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repeated.  If  no,  the  number  of  survivors  and  number  of  air- 
craft shot  dovm  are  printed  and  the  simulation  ended  by  set- 
ting MSTOP  to  a negative  value. 


If  the  aircraft  passes  through  the  threat  area  without 
being  fired  on  or  if  the  threat  does  not  fire  all  of  its  mis- 
siles, then  subroutine  AORNG  is  called  when  the  aircraft 
crosses  MINEFR.  A listing  of  this  subroutine  followsi 


ATRIB(1)=1.E20 
ATRIBU)  = 2.0 
CALL  FILEM(l) 
ATRIB(1)=1.E20 
ATRIB(2)=3.0 
CALL  FILEM(l) 
CALL  OTPUT 
110  RETURN 


This  subroutine  stops  SCAN  and  SCANT  from  searching  for  the 
aircraft  by  setting  their  time  of  next  occurrence  far  into 
the  future.  A call  is  then  made  to  OTPUT. 


This  is  the  last  of  the  Fly  By  programs.  The  results 
of  FB4  are  discussed  in  Chapter  V.  aind  listed  in  Appendix 


Each  missile  system  has  an  engagement  envelope.  This 
envelope,  ideally,  encompasses  the  area  around  the  site 
where  the  missile  could  be  used  to  kill  an  aircraft.  A 
general  picture  of  this  envelope  is  shown  in  Fig.  9 below. 


RANGE  ALTITUDE 

Fig.  9.  Ideal  Missile  Engagement  Envelope 


The  missile's  characteristics,  combined  with  those  of  the 
tracking  and  guidance  device,  determine  its  maximiim  slant 
range,  minim\jm  slant  range,  and  maximum  and  minimum  altitude 
of  effectiveness.  These  parameters  give  us  the  "ideal"  mis- 
sile envelope.  But  there  is  more  to  the  problem. 

A slow  missile  would  not  catch  an  aircraft  going  away. 
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Also,  the  tracking  device  requires  a certain  amount  of  time 
to  determine  a firing  position,  even  if  the  aircraft  is  al- 
ready within  range.  Infra-red,  or  heat-seeking,  missiles 
are  less  effective  when  used  against  on-coming  aircraft. 
These  parameters  and  their  limitations  lead  to  an  operation- 
al envelope.  The  missile  is  fired  only  when  there  is  a rea- 
sonable chance  of  success.  A general  picture  of  what  an 
operational  envelope  might  look  like  is  shown  in  Fig.  10 
below. 

range  altitude 


max  ceiling 


(0,0,0)  (0,0,0) 


Fig.  10.  Operational  Engagement  Envelope 

Most  of  the  parameters  described  above  have  been  used 
in  the  Fly  By  programs.  Some  of  these  are  missile  speed, 
XMSPED,  its  minimum  and  maximum  effective  ranges,  MINERF 
and  EFFRM,  and  the  effective  range  of  the  tracking  device, 
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EFFRT.  The  numbers  used  to  denote  these  parameters  in  the 
Fly  By  programs  are  not  peculiar  to  any  one  system.  To  keep 
the  report  unclassified,  they  were  simply  selected  arbitrar- 
ily but  kept  reasonable.  The  program  is  flexible  so  that 
each  of  these  parameters  can  be  changed  to  reflect  specific 
systems  in  the  future.  It  would  also  be  rather  easy  to  add 
other  parameters  necessary  to  make  the  programs  more  realis- 
tic. 

One  other  device  of  interest  in  this  engagement  is  the 
optical  device  being  used  by  the  SAM  guidance  or  tracking 
system.  Data  on  exactly  what  happens  to  optical  systems 
when  hit  by  various  known  or  proposed  countermeasures  de- 
vices is  currently  being  sought.  Some  experiments  are  al- 
ready being  conducted  and  others  are  in  the  planning  stages. 
These  optical  devices  have  not  been  incorporated  into  the 
Fly  By  programs.  However,  when  the  data  has  been  obtained 
and  the  pertinent  parameters  determined.  Fly  By  4 will  re- 
quire few  or  no  changes  to  incorporate  the  device  as  the 
tracker  in  the  SCANT  and  TRACK  subroutines.  What  data  that 
is  obtainable  is  being  gathered  at  the  AFAL  library,  where 
it  will  be  available  for  future  use. 

Other  parameters  of  the  tracking  device  need  to  be  in- 
corporated to  make  the  program  realistic.  For  an  optical 
device,  the  contrast  of  the  aircraft  plays  an  important  role 
The  characteristics  of  the  atmosphere  will  also  have  to  be 
taken  into  account.  The  probability  of  lock-on  must  be  made 
to  reflect  these  parameters. 
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OCM  Parameters 
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A SAM  will  kill  if  it  reaches  a point  in  the  air  where 
the  warhead  is  within  lethal  range  of  an  aircraft  and  is  then 
detonated.  This  kill  distance  varies  with  the  warhead  type. 
The  essence  of  a countermeasure  device  is  to  degrade  the 
tracking  device  to  the  point  that  it  cannot  guide  a missile 
within  lethal  range  of  the  aircraft.  This  could  occur  after 
a missile  has  been  launched,  since  most  missiles  are  guided 
to  the  target.  For  AAA  systems,  this  is  not  the  case.  Once 
the  shell  leaves  the  gun,  it  does  not  matter  that  the  track- 
er cannot  see  you,  at  least  not  as  far  as  that  shell  is  con- 
cerned. So  in  the  case  of  optically  guided  AAA,  you  must 
stop  the  tracker  from  ever  knowing  your  position  accurately 
enough  for  the  AAA  to  become  effective.  Thus,  the  purpose 
of  the  countermeasures  device  is  to  introduce  errors  into 
the  tracking  system  of  the  threat  and  to  thus  degrade  the 
system  to  such  a point  that  its  effectiveness  is  greatly 
reduced. 

No  data  was  available  on  the  devices  under  current  study 
to  be  used  as  optical  countermeasures  systems.  The  param- 
eters of  these  devices  are  known,  but  are  sensitive  and 
cannot  be  discussed  in  this  report. 


Fly  By  ^ 

The  next  step  in  developing  the  optical  analysis  is  the 
Fly  By  5 program.  This  program  will  incorporate  new  sub- 
routines and  changes  to  old  ones  to  reflect  a realistic 
optical  engagement.  A discussion  of  changes  will  be  first. 
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In  the  Fly  By  programs  in  this  report,  the  detection 


and  tracking  probabilities  are  generated  using  the  range  of 
the  aircraft  and  the  number  HIALT.  This  number  incorporates 
all  the  parameters  not  specifically  mentioned,  such  as  air- 
craft cross-section,  contrast,  power  of  the  radar  or  optical 
device,  and  field  of  view  of  the  tracker,  into  one  n\Amber 
which  is  then  used  to  determine  the  probabilities.  This  is 
too  simplified  for  a realistic  engagement.  The  HIALT  n\imber 
must  be  calculated  using  the  above  parameters  and  any  others 
which  might  apply.  There  are  several  computer  programs  in 
the  RF  regime  which  could  be  used  as  guides  in  changing  the 
INTLC  subroutine  to  reflect  the  new  HIALT.  One  of  these  is 
the  pool  AAA  program  used  at  AFAL.  It  was  used  extensively 
to  obtain  an  overall  view  of  the  engagement  process  described 
in  this  report  (Ref  2) . Another  approach  would  be  to  write 
a new  subroutine  called  PROBD  to  determine  the  value  of  the 
number  HIALT.  Either  method  would  give  a more  realistic 
tracking  and  detection  model. 

Other  changes  would  have  to  be  made  in  the  FI REM  sub- 
routine to  reflect  the  possible  use  of  different  detonating 
devices.  A missile  being  detonated  by  an  observer  on  the 
ground  or  at  a computer  calculated  time  of  closest  approach 
would  have  a different  effect  than  one  with  a contact  or 
proximity  fuse.  Also,  some  missiles  can  be  fired  at  aircraft 
going  away.  This  would  have  to  be  incorporated  into  the 
FIREM  and  TRACK  subroutines. 

To  reflect  countermeasure  degradation  of  the  tracking 
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system,  a new  subroutine  JAMMER  must  be  written.  As  envi- 
sioned, this  subroutine  would  incorporate  the  parameters  of 
the  OCM  devices  and  deliver  as  output  a tracking  degradation 
number.  This  number  would  be  used  in  the  TRACK  subroutine 
to  determine  where  the  tracker  sees  the  aircraft.  All  of 
the  parameters  of  the  OCM  device  would  be  involved,  as  well 
as  the  atmospheric  properties.  Input  would  come  from  another 
new  subroutine,  RECEIV. 

RECEIV  would  model  the  countermeasures  receiver  and  its 
ability  to  detect  optical  threats.  The  parameters  of  the 
threat,  the  aircraft,  and  the  atmosphere  would  be  combined 
to  deliver  the  range,  azimuth,  and  priority  of  the  threat. 
This  data  would  be  sent  to  JAMMER.  A flowchart  showing  where 
the  two  new  subroutines  would  be  incorporated  into  Fly  By  4 
is  shown  in  Fig.  11  on  page  ^0, 
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FROM  FLY 

BY  4, 

PAGE  33. 

TTIME 

>10.0? 

YES 

TO 

FLY  BY 

4 

34 

PAGES 

32 

/ TRACK 


LOCK-ON? 


' IS  THE  > 
ACFT  WITHIN 
. EFFRM?  > 


CALL 

RECEIV 


/WAS  \ 
/ LOCK-ON  ^ 
DETECTED  BY 
\ ACFT?  y 


COMPUTE  INTERCEPT 
POINT 


DETERMINE  POSITION 
OF  ACFT  AS  SEEN  BY 
TRACKER  


CALL  JAMMERM — ^DEGRADE  TRACKII 


Fig.  11.  Flowchart  of  FB4  Changes  for  FB5 


V.  Results 

General 

Each  Fly  By  program  has  a different  set  of  data  cards. 
To  help  keep  track  of  the  values  used  by  GASP,  each  program 
prints  out  a simulation  project  number  followed  by  all  the 
values  read  in  on  the  GASP  data  cards.  A detailed  explana- 
tion of  these  values  can  be  found  in  Ref  1,  Chapter  3. 

Another  GASP  printout  lists  the  contents  of  all  the 
files  and  the  contents  of  the  state  variable  storage  area. 
This  is  done  twice  for  each  program.  The  headings  are  GASP 
FILE  STORAGE  AREA  DUMP  and  GASP  STATE  STORAGE  AREA  DUMP. 

The  value  of  TNOW  at  the  time  of  the  dump  is  listed  after 
the  title  in  the  printout.  A dump  is  made  for  values  of 
TNOW  at  the  beginning  and  end  of  the  simulation. 

One  other  GASP  printout  occurs.  This  is  the  GASP  sum- 
mary report,  which  lists  the  project  nxmiber,  date,  parameter 
sets,  and  the  statistics  collected  on  each  variable.  All  of 
these  printouts  are  listed  with  the  user  results  in  Appendix 
A. 

Fly  By  1 

Project  1 is  the  simplest.  The  aircraft  is  detected 
as  it  crosses  TEFFR  inbound  and  outbound.  No  statistics  sure 
collected  and  only  one  aircraft  is  involved.  The  INTERMED- 
IATE RESULTS  are  as  follows i 

TNOW  = 6.00 

COORDINATES  ARE 
RANGE=  99^7 . 86 


L 


0.00  8600.00  5000.00 


51 


TNOW=  49 . 25 

COORDINATES  ARE 
RANGE=  10034.44 


0.00  -8700.00  5000.00 


i 

j 


1 

1 


As  can  be  seen,  the  aircraft  crossed  TEFFR  inboxind  at  time  j 

I 

TNOW=6.00s.  The  actual  range  was  9947.86m,  which  is  well  j 

i 

within  the  200m  tolerance  set  in  KROSS.  As  the  aircraft 
left  the  area,  it  crossed  TEFFR  at  TN0W=49.25s  at  a rsinge 
of  10034.44m,  again  well  within  tolerances. 

Althoxigh  the  results  are  rather  simplistic,  one  of  the 
major  purposes  of  Fly  By  1 was  to  become  familiar  with  the 
GASP  executive.  Fly  By  1 was  accomplished  with  this  in  mind 
and  succeeded  in  laying  the  groundwork  for  continuing  to 
Fly  By  2. 

Fly  ^ 2 

In  Fly  By  2,  the  SCAN  and  DETECT  subroutines  were  added 
to  make  the  detection  process  more  realistic.  Again,  no  sta- 
tistics were  collected.  The  INTERMEDIATE  RESULTS  are  listed 
be low t 

THE  ACFT  IS  IN  RANGE 
TNOW  = 6.00 

COORDINATES  ARE  0.00  86OO.OO  5OOO.OO 

RANGE=  9947 . 86 

THE  ACFT  HAS  BEEN  DETECTED 

TNOW  6.00 

COORDINATES  ARE  0.00  86OO.OO  5OOO.OO 

RANGE=  9947.86 

The  flags  about  range  and  detection  indicate  that  the  infor- 
mation following  concerns  time,  coordinates,  and  range  of 
crossing  TEFFR  and  detection  by  SCAN,  respectively.  The 
values  used  for  the  different  variables  caused  detection  to 
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occur  when  the  aircraft  crossed  TEFFR.  The  information  print- 
ed in  Fly  By  1 about  aircraft  withdrawal  does  not  appear  be- 
cause the  simulation  ends  when  the  aircraft  is  detected.  If  ’ 

the  aircraft  was  not  detected,  the  AORNG  subroutine  would 
print  out  the  aircraft  withdrawal  information  beneath  the 
flag  "THE  ACFT  IS  OUT  OF  RANGE".  A complete  listing  of  the 
Fly  By  2 results  and  printout  begins  on  page  63  in  Appendix 
A. 

IllMl 

This  program  is  the  first  in  which  more  than  one  air- 
craft was  used  and  statistics  collected.  Use  of  more  them 
one  aircraft  does  not  mean  the  one-on-one  has  been  abandoned, 
but  that  it  is  being  repeated  several  times.  This  is  re- 
flected throughout  the  results  and  printout,  which  is  listed 
in  Appendix  A beginning  with  page  66.  The  first  few  lines 
of  the  INTERMEDIATE  P?ESULTS  are  listed  belowi 

THE  AIRCRAFT  HAS  BEEN  DETECTED 
TNOW  = 73.75 

COORDINATES  ARE  0.00  27250.00  5OOO.OO 

RANGE  IS  27704.92 

THE  AIRCRAFT  HAS  BEEN  DETECTED 
TNOW  = 15.75 

COORDINATES  ARE  0.00  3885O.OO  5OOO.OO 

RANGE  IS  39170.43 

All  twenty  aircraft  are  detected  before  passing  out  of  the 
threat  area  and  the  appropriate  flag,  time,  position,  and 
range  of  each  detection  are  listed.  The  data  above  repre- 
sents two  aircraft.  In  addition,  statistics  are  collected 
and  printed  in  the  summary  report.  A listing  of  these 
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statistics  is  shovm  below 


RNG  DET 


MEAN 

STD  DEV 

SD  OF  THE  MEAN 

CV 

MINIMUM 

MAXIMUM 

OBS 


.369IE+O5 

.3092E+O4 

.6915E+O3 

.8377E-01 

.277OE+O5 

.3996E+O5 

20 


The  closest  an  aircraft  came  to  the  site  before  detection 
occurred  was  27250.00in.  The  mean  range  of  detection  was 
36910.0m,  with  a standard  deviation  of  691.5® 

The  last  data  printed  out  by  Fly  By  3 is  a histogram 
of  the  range  of  detection.  It  is  shown  in  Appendix  A on 
page  72.  This  histogram  shows  relative  and  cumulative  fre- 
quency of  observation  of  the  values  between  the  upper  and 
lower  limits  of  each  cell.  The  C indicators  are  for  the 
ciomulative  graph,  while  the  * indicators  tell  what  percent- 
age of  the  total  observations  fell  within  each  range  of 
values. 


Fly  By  4 

In  Fly  By  4,  the  flags  for  aircraft  detection,  penetra- 
tion, and  withdrawal  were  dropped.  Therefore,  the  INTER- 
MEDIATE RESULTS  only  contains  information  on  aircraft  sur- 
vivability. These  results  are  listed  belowi 

THE  NUMBER  OF  SURVIVING  ACPT  IS  6 
THE  NUMBER  SHOT  DOWN  IS  l4 


Data  has  also  been  collected  on  three  variables.  The 
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statistics  are  calculated  and  then  printed  out  in  tabular 
form.  They  are  listed  belowi 


mAH 
STD  DEV 
SD  OF  MEAN 
CV 

MINIMUM 

MAXIMUM 

OBS 


RNG  DET 

.1422E+06 

.1364E+04 

.305IE+O3 

.9593E-02 

.1379E+06 

.1432E+06 

20 


RNG  TRK 

.5024E+05 

.2948E+04 

.6592E+O3 

.5868E-OI 

0.0 

.538OE+O5 

20 


RNG  FIR 

.3952E+05 

.5463E+03 

.951OE+02 

.1382E-01 

0.0 

.3995E+O5 

33 


The  printout  is  completed  with  the  histograms  of  the  three 
variables.  These  are  listed  in  Appendix  B,  starting  on 
page  . As  would  be  expected,  a large  number  of  the  air- 
craft are  detected  as  soon  as  they  enter  TEFFR.  the  RNG  TRK 
histogram  shows  that  the  tracker  "locked-on"  to  all  but  one 
aircraft  as  they  came  within  range.  And  once  "locked-on", 
the  threat  fired  a missile  at  each  aircraft  as  it  crossed 
EFFRM.  Seven  aircraft  were  destroyed  by  the  missiles  on 
the  first  shot,  while  only  six  survived  the  second  round. 

Summary 

With  Fly  By  4,  a point  has  been  reached  at  which  a true 
optical  engagement  can  now  be  realistically  programmed.  The 
program  scans,  detects,  tracks,  and  fires  at  the  aircraft 
with  parameters  that  are  set  by  the  user.  With  this  flexi- 
bility and  the  addition  of  subroutines  representing  the  re- 
ceiver and  jammer  on  the  aircraft.  Fly  By  5 will  give  some 
realistic  data  on  the  effectiveness  of  optics.  The  limita- 
tions of  the  program  have  been  noted,  but  the  flexibility 
allows  the  program  to  easily  model  any  situation. 
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VI.  Recommendations 


The  ultimate  goal  of  this  project  is  to  build  a computer 
program  to  analyze  a one-on-one  engagement  in  the  optical 
regime.  As  stated  before,  certain  additions  to  the  Fly  By 
4 program  must  be  made  to  reach  this  goal.  Data  on  the  re- 
ceiver and  countermeasures  devices  must  be  collected  and 
analyzed  to  determine  what  parameters  need  to  be  included  in 
the  JAMMER  and  RECEIV  subroutines.  This  means  a large  amount 
of  data  must  be  collected.  One  possible  way  would  be  to 
propose  an  AFIT  thesis  to  research  and  collect  data  on  opti- 
cal receivers  and  jammers  and  put  it  into  usable  form.  A 
tabulation  of  these  devices,  both  existing  and  proposed, 
would  be  invaluable  to  individuals  in  the  AFAL  who  will  con- 
tinue this  project.  Many  of  the  docxaments  necessary  to  com- 
plete such  an  analysis  have  been  ordered  by  the  AFAL  library 
but  were  not  available  when  this  report  was  written.  It  is 
felt  that  collection  and  compilation  of  data  is  essential 
to  the  successful  completion  of  the  one-on-one  project. 

Another  area  requires  the  collection  of  large  amoimts 
of  data.  The  characteristics  of  the  threat  optical  devices 
need  to  be  researched  and  cataloged.  Although  many  sources 
of  this  data  exist,  only  a few  of  them  are  at  the  AFAL  li- 
brary. Those  ordered  for  this  report  were  not  available  at 
its  completion.  However,  they  could  not  have  been  included 
if  this  report  was  to  remain  unclassified.  The  data  exist 
and  must  be  collected  and  tabulated. 

One  large  grey  area  which  was  discovered  during  this 
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report  concerns  the  effect  of  optical  countermeasures  devices 
on  optical  threat  devices.  Several  experiments  have  been 
conducted  or  are  being  conducted  at  this  time  concerning 
these  effects.  Information  available  indicates  wide  var- 
iances in  the  results.  The  pertinent  questions  include  the 
followingj 


1.  What  effect  does  the  countermeasures 
device  have  on  the  optics  and  its  ability 
to  see  the  aircraft? 

2.  How  much  error  must  be  introduced  to  give 
reasonable  chance  of  aircraft  survival? 

3.  How  much  bloom  time  is  required  to  intro- 
duce a specific  error? 

4.  Does  a laser  device  destroy  a video  tube 
or  just  cause  it  to  bloom? 

5.  What  power  level  is  required  for  destruc- 
tion? 

6.  What  power  level  is  required  for  blooming? 

7.  How  long  does  blooming  last  at  different 
power  levels? 


These  and  other  questions  must  be  answered  and  put  into  a 
form  which  can  be  included  in  the  parameters  of  Fly  By  5* 

An  AFIT  thesis  on  the  effects  of  a laser  hitting  a TV  video 
tube  would  go  a long  way  toward  answering  these  questions 
and  would  give  AFAL  a key  on  which  to  evaluate  other  data. 
This  project  will  be  proposed  to  the  current  AFIT  physics 
class. 
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Information  marked  by  * is  GASP  generated.  A detailed 
explanation  can  be  found  in  Ref,  1. 
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Appendix  B 

Calculation  of  the  Intercept  Point 

Variables 

t = missile  flight  time 

V = velocity  of  the  aircraft 
tJ  = velocity  of  the  missile 

V = speed  of  the  aircraft 
u = speed  of  the  missile 

Rjj  = initial  position  of  the  aircraft 

r^  = initial  position  of  the  missile 

R,r  = position  of  the  aircraft,  missile  at  intercept 


Diagrams 
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Solution 


The  initial  positions  of  the  aircraft  and  missile  are 
given  by 

^0  = *0^  * ^0^  ■"  ®o’' 

In  a given  time  t,  the  aircraft  and  missile  will  move 
to  new  positions  defined  by 

R = Rq  + 7t 

r = ? + Ut 

o 

For  intercept  to  occvir, 

R = r 

Therefore, 

Ro  ^ = ^o  ^ Ut 

or 

Ro  - ^o  = 

Breaking  these  into  components,  we  have 

Xi  + Yj  + Zk  - Xpi  - y^jj  - Zpk  = 7t  - iTt 

j 

P 

But  the  velocities  can  also  be  broken  up  into  components, 

I 

: 80 

I 

I 

I 
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Vt  - Ut  = (V^i  + Vyj  + V^k  - - Uyj  - U^k)t 


(Xq  - Xq)(-u)(1  - cos^e)^/^  = (V  - ucos0)(Zq  - Zq) 


Squaring, 


^^0  *0^  -(u^)(l  - cos^e)  = - 2VUCOS9  + u^cos^e 


(2o  - 


(X.  - x„) 


o ■ ,.2 


(Zq  - 


(u")  - 


^^0  ■ *0^  /.  2_  2 


(Zo  - 


(u  cos  9)  = V - 2VUCOS9 


2 2 
+ U^C0S^9 


((  + u^)(cos^9)  - 2VUOOS9  + 


(Z„  - 


“o)' 


<^0  - "o)' 


Solve  for  cos9  using  the  quadratic  formula,  with  the 
following  constants 

-b  = 2Vu 


b'”  = 


a = 


c = 


(Zq  - J^q)' 
'^0  - *o)' 


rU  + U 


- 


(^0  - 
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